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Abstract: With the advent of Talbot-Lau interferometers for x-ray phase- 
contrast imaging, oblique and grazing incidence configurations are now 
used in the pursuit of sub-micron grating periods and high sensitivity. Here 
we address the question whether interferometers having oblique incident 
beams behave in the same way as the well-understood normal incidence 
ones, particularly when the grating planes are non-parallel. We derive the 
normal incidence equivalence of oblique incidence geometries from wave 
propagation modeling. Based on the theory, we propose a practical method 
to correct for non-parallelism of the grating planes, and demonstrate its 
effectiveness with a polychromatic hard x-ray reflective interferometer. 
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1. Introduction 

Lau and Talbot-Lau grating interferometers produce wave interference effects from 
polychromatic, spatially incoherent sources and thus have wide ranging applications in 
molecular, atomic [1,2] and x-ray interferometry [3-6]. A Talbot-Lau interferometer consists 
of two planar gratings and a detector plane which are arranged in parallel (Fig. 1) [7-9]. 
These are the three principal planes of the interferometer. When certain symmetry conditions 
are met the interference patterns from different wavelengths and different incident angles are 
phase-locked at the detector plane, which enables interferometry with low coherence, 
polychromatic sources [1,2,8,10]. 

Although the literature on Talbot-Lau grating interferometers have dealt with normal 
incidence geometry, recent efforts in developing sub -micron grating periods for high- 
sensitivity x-ray interferometers have led to oblique incidence gratings, such as multilayer- 
coated transmission gratings [11] and grazing-angle reflective gratings [12-17]. Now there is 
a need for a theoretical description of oblique incidence configurations under incoherent 
illumination. Of particular interest is the effect of slight non-parallelism of the three principal 
planes, which becomes significant with oblique incidence. Here we will give a first-principle 
derivation of the equivalent normal incidence geometry for oblique incidence geometry under 



Detector plane 



Light source 




Fig. 1. An illustration of a Talbot-Lau grating interferometer. The two gratings and the detector 
plane are parallel. The light source can be polychromatic with low spatial coherence and still 
produces interference fringes on the detector plane if the appropriate grating periods and 
geometric layout are used. 
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either parallel beam or cone beam illuminations, and obtain the leading order effects of non- 
parallelism of the three principal planes of the interferometer. Point sources will be 
considered first and the results will be applied to incoherent area sources. We will show that 
the effect of slight non-parallelism is equivalent to changing the inter-plane distances in a 
parallel configuration. The practical implications of the above understanding on oblique- 
incidence transmission and grazing angle reflective interferometers will be discussed. We 
demonstrate with a hard x-ray reflective interferometer that non-parallel grating planes cause 
a loss of interference fringe visibility, which can be recovered with distance adjustments. 

2. Basic theory 

The basis of the derivation is the Fresnel-Kirchhoff principle [18], which states that the wave 
function at a point on the detector plane arises from the superposition of waves proceeding 
from all points on the exit surface of the Gl grating, and the wave function at a point on the 
entry surface of the Gl grating arises from the superposition of waves proceeding from all 
points on the exit surface of the GO grating, and so forth. Figure 2 shows a generalized layout 
of the interferometer. For the sake of clarity we confine the derivation to the case where the 
beam axis and the normal vectors of the three principal planes are coplanar, and call this plane 
the pitch plane. 




Fig. 2. Coordinate system and general interferometer configuration described by the theoretical 
derivations. The pitch angles of the 3 principal planes relative to the axis of the incident beam 
are 0, 8j and 0 2 . The coordinates within the principal planes are (x„, y n ) and the Y axes of the 
planes are all perpendicular to the incident beam axis and parallel to each other. The para-axial 
distances between the source and the 3 principal planes are R,. In the case of a parallel incident 
beam the distance R 0 is set to infinity. 

In the general case illustrated in Fig. 2, the incident angle 8 onto the GO grating is defined 
as the angle between the central axis of the cone beam and the GO plane, and so on for the 
other planes. Several configurations will be looked at in order of increasing complexity, 
including CI the simplest case of a parallel beam illuminating three parallel planes at normal 
incident angle; C2 a cone beam from a point source illuminating three planes on axis at 
normal incident angle; C3 a cone beam illuminating parallel planes off-axis at an oblique 
incident angle; C4 the most general configuration of a cone beam illuminating 3 non-parallel 
planes at an oblique incident angle. We will show the equivalent CI configuration of a given 
C2 configuration, then the equivalent C2 of a given C3, and lastly the equivalent C3 of a 
given C4 when the differences among the pitch angles of the three planes are small. 
Following this procedure, all configurations are mapped to the simplest parallel beam normal 
incidence configuration. For clarity, the symbol "~" will be used to relate two functions that 
differ by a constant multiplication factor, and such factors are left out when they do not affect 
the essential results of the derivation. 
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The general C4 configuration illustrated in Fig. 2 is defined by a set of inter-plane 
distances R h incident angles 0„ and grating transmission functions G m . The Kirchhoff integral 
of the wave amplitude A 4 on a point (x 2 , y 2 ) on the detector plane is 



A 4 (x 2 ,y 2 \Rj,e n ,G m )~ j" dx 0 dy 0 dx x dy, 



exp(-ife|r () |) 
— u 

l r o| 



(l) 



exp(-/fc|r,|) 



G l (x l ,y 1 )(n l -e 2 



) 



exp(-/A: r 2 ) 



r, 



where k is 27i/wavelength, (x 0 , y 0 ) and (x y , y } ) are coordinates on the GO and Gl grating planes 
as illustrated in Fig. 2, r 0 , r : and r 2 are vectors from the source point to (x 0 , yd), from (x 0 , y 0 ) to 
(x y , yj), and from (x y , yj) to (x 2 , y 2 ), respectively, ei and e 2 are the directional unit vectors of r : 
and r 2 , G 0 (x 0 , y 0 ) and Gy(x y , y^) are the transmission functions of the GO and Gl gratings, 
respectively, and n 0 and ni are the normal vectors of the GO and Gl planes, respectively. 

To proceed further the Fresnel approximation will be used, which is valid under the 
assumptions that the diffraction angles from the GO and Gl gratings are small. The 
approximations are i) The inclination factors (n 0 -ei) and (n^e^ in Eq. (1) are constants 
independent of the coordinates (x 0 , y 0 ) and (x h y y ), and the denominators in Eq. (1) are 
approximated to the para-axial distances between the source and the planes, R 0 , R 1 and R 2 , as 
illustrated in Fig. 2. ii) The lengths of the r 0 , r l and r 2 vectors in the phase delays are 
approximated to 2nd order terms of the coordinates in the grating and detector planes, iii) The 
effects of non-parallelism of the planes are considered only to the leading order terms of the 
differences between the incident angles on the planes. 

Under the above approximations, the wave amplitude of the C4 configuration is simplified 
from Eq. (1) to 



The key quantity that varies among the different interferometer configurations is the phase 
factor in the above expression, which is proportional to the sum of the propagation distances 
between the planes. To show that one configuration is equivalent to another, it is necessary to 
show that the phase factors of the two configurations differ by only a constant for the specific 
pair of the layout parameters (Rj, 6>„, G m ) that describe the two configurations. To do so, the 
phase factor of the general C4 configuration will be expanded to 2nd order terms of the 
coordinates in the 3 planes, then simplified expressions will be derived for the less general CI 
to C3 configurations. It will then be possible to find the pair of layouts (Rj, Q n G m ) for the C2 
and CI configurations that result in equivalent expressions of the phase factors, and do the 
same between C3 and C2 and so forth. By this procedure, all configurations are ultimately 
mapped back to the simplest parallel beam normal incidence interferometer. 

Denoting the differences between the incident angles 8 h 8 2 relative to 6 as 



then under the approximations ii and iii, the lengths of the vectors r 0 , r ; and r 2 in Eq. (2) are 
expanded to the second order terms in the X and Y coordinates and first order terms in the S 
angles: 



A 4 (x 2 , y 2 1 Rj , 6 n , G m ) ~ j dx 0 dy 0 dx l dy 1 



1 



" u 1 1 p j? i? 
expH'fc(|r 0 1 + |r, | + |r 2 |)]G 0 (x 0 , y 0 )G, (x, , y, ) 



(2) 



|r 0 1 = [(R 0 cos 6 + x () ) 2 + (R 0 sin 6) 2 + y 0 2 ] 




(4) 
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|r, | * [(/?, cos 9 + x, - x 0 ) 2 + (R, sin 9 - x^f + (y, - y Q ) 2 ] 1/2 

« *, + (x, - * 0 ) cos 0 - M sin 0 + i (Xl - X " )2sin ^ + I (y '- y " )2 (5) 



(x, -x^x^, sin#cos6> 



|r 2 | » [(/? 2 cos# + x 2 -Xj) 2 + (/? 2 sin6>-x 2 5 2 +x l( ?,) 2 +(y 2 - y^ 2 ]" 2 

« /? 2 +(x 2 -x,)cos6>-(x 2 <? 2 -Xj(?,)sin^ (6) 

1 (x 2 -x,) 2 sin 2 9 1 (y 2 - y^ 2 (x 2 -Xj)(x 2 5 2 -x^sin^cosf? 

2 R 2 2 R 2 R 2 

The phase factor in Eq. (2) is proportional to the sum of the above 3 lengths. Summing up 
the expressions in Eqs. (4)-(6) yields an expression for the phase factor of the C4 
configuration: 

2 2 

<f> 4 (x 2 , y 2 I Rj , 9 n , GJ = -*£ |r, | ~ £ Rj + (cos 9 - S 2 sin 9)x 2 

0 0 

1 . 2 /ir 2 s 1 1 \ 2/1 1 \ -^7 2x n X, 2X.X-, 

+ - sin 6>[x 0 ( — + — ) + x, ( — + — ) + — — — ] 

2 R Q R l R l R 2 R 2 R l R 2 



. . „ 2 ,,1 1, x 2 8 2 x 0 x 1 S l x.xJS, + S 2 ) 1 
+ sin 9 cos 9[x 1 S l ( — + — ) + - 2 — — — L ^ J — ] 



(7) 



R i R 2 R 2 R i R : 



1 r 2, i , K , i, i , i , , y 2 2 y 0 yi 2 yj 



+-[y 0 (—+—) + y, (—+—)+ ., , 

2 /?j /? 2 /? 2 /?j /? 2 

Since in the expression of Eq. (1) the integrations are over the (x 0 , y 0 ) and (x lt y t ) 
coordinates, the terms in the phase factor that are independent of the two coordinates will be 
lumped into a quantity tp c (x 2 , y 2 ), and Eq. (7) can be written as 

<f> 4 (x 2 ,y 2 \R p 9 n ,G m )~<l> c {x 2 ,y 2 ) 

1 . 2 rw 2 / 1 1 \ 1\ 2x n x, 2x.x-, 
+ -sin 6>[x 0 ( — + — ) + xf ( — + — ) — — ] 

2 R 0 R^ R^ R 2 R l R 2 

n nr 2? A K X 0 X A X 1 X 2 (^1 + <^2 ) n ^ 

+ sin 6> cos 9[x l S l ( — + — ) — 9 - LJ - — LJ — ^ — ] 

R^ R 2 R\ R 2 

4[ yo 2 (- + -) + ^(- + -)-^-^], 

2 /?„ R { R 2 /? 2 



where 



<j> c (x 2 , y 2 ) = ^ + (cos # - J 2 sin #)x 2 



+ I Sm 2^ +Sm 0 COS 0 X 2 <5 2 , 1 ^2 



(9) 



R 2 R 2 2 R 2 



Generally the detector only detects the intensity of the wave, and so the phase term tp c (x 2 , 
y 2 ) has no effect on the detected pattern and can be neglected in further derivations. 

Additionally, the tilt of the principal planes are confined in the XZ plane (Fig. 2), and 
therefore in the phase factor expression of Eq. (8), the terms involving the Y coordinates are 
not affected by the oblique incident angles. Therefore in searching for the equivalent 
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configurations we will first attempt to equate the terms involving the X coordinates, and then 
apply trivial scaling to the Y coordinates accordingly. 

3. From parallel beam normal incidence to cone beam normal incidence 

These two configurations (CI and C2) are defined by 

CI : (tf 0 =«>,0„ = 90°, GJ, (10) 

and 

C2: (Finite R,,0 n = 90°, GJ. (11) 

The phase factor of the wave on the detector plane, expressed in general in Eq. (8), is now 
simplified to 

A (x 2 , y 2 I R 0 = oo, 6> = 90°, G m ) ~ <f> c (x 2 ,y 2 ) 

+ -[x 0 2 - + x 1 2 (i + -)-^-^] (12) 
2 R { R 2 R^ R 2 

+ -[y 0 — + yi (— + — ) — — 
2 R x R 2 R t R 2 



and 

</> 2 {x 2 ,y 2 \R p G n = 90°,GJ~^(x 2 ,y 2 ) 

+ -[x () 2 (— + -) + x 1 2 (- + — )-^-^] (13) 
2 /?q /?j R 2 Ry R 2 

+ -u 2 (- + -) + y 1 2 (- + -)-^-^]. 

2 Rq R^ R^ R ? R^ R 2 

To find the equivalent CI configuration for C2, the terms that contain the x 0 and x y 
coordinates in Eq. (13) should equal the corresponding terms in Eq. (12). Denoting the 
parameters in the equivalent CI configuration with the above condition means that 

x' 2 11 

^ = x () 2 (— + —), (14) 
K i Kq K { 

(15) 
(16) 



r: r, 



R 2 R 2 



1 - ,2, 1 1 



jc'i (— + —) = x'i (— + —), (17) 

1 /?; /? 2 1 /?, /? 2 

and the same set of equations apply to the Y coordinates. Although the solutions to these 
equations are by no means unique, one can limit them to a single solution by requiring that the 
Gl grating transmission function is not altered. In other words, 

x\ = x r (18) 
Substituting Eq. (18) into Eq. (15) yields 

^JL = ^L. (19) 
R\ R, 
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Combining Eq. (19) and Eq. (14) yields the rescaling relationships of 

R'^R^l + R./R,), (20) 

and 

jc' 0 =x 0 (l + R l /R 0 ). (21) 
Substituting Eq. (18) and (20) into Eq. (17) yields the scaling relationship forR 2 : 

R' 2 =R 2 R ° +Rl . (22) 
Rq + R^ + R 2 

Substituting Eq. (18) and (22) into Eq. (16) yields the scaling for the x 2 coordinate 

R " 4 (23) 



Rq + 7?j + 

Thus, in the equivalent parallel beam configuration to a cone beam geometry, the spacing 
between the principal planes are rescaled according to Eq. (20) and (22). The rescaling of the 
X coordinates given in Eq. (18), (21) and (23) also applies to the Y coordinate. 

The grating transmission functions are rescaled according to the general relationship 

G' m (x' m ,y' m ) = G m (x m ,y m ). (24) 

Based on Eq. (21) and (24), the GO grating transmission function of the equivalent parallel 
beam configuration is related to the cone beam one by 

G ' 0 (x 0 , y 0 ) = G 0 (x 0 /(1 + RJ R 0 ), y 0 /(l + RJ R 0 )). (25) 

By choice the Gl grating transmission function is unchanged: 

G\(x l ,y 1 ) = G l (x l ,y 1 ). (26) 

Lastly, the patterns on the detector planes of the two equivalent configurations are related 
to each other by 

I 2 (x 2 ,y 2 )~I\(x\,y\), (27) 

where I 2 is the detected intensity pattern of the C2 configuration, and Ij ' the detected pattern 
of the equivalent CI configuration. Based on Eq. (23), this relationship is explicitly 

I 2 (x 2 ,y 2 )~I\( R ° +Rl x 2 , R ° +Rl y 2 ). (28) 
R^ + R^ + R 2 Rq + R { + R 2 

In summary, the relationships expressed in Eqs. (20), (22), (25), (26) and (28) completely 
describes an equivalent parallel beam, normal incidence interferometer for a cone beam, 
normal incidence interferometer as long as the Fresnel approximation holds. 

4. From cone beam normal incidence to cone beam oblique incidence 

The next step involves the cone beam oblique incidence (C3) configuration, which is defined 

as 

C3 : (finite R„G n =6,G m ). (29) 
For C3 the phase factor of the wave on the detector plane is simplified from Eq. (8) to 

</> i (x 2 ,y 2 \R j ,0 n =0,GJ~</> c (x 2 ,y 2 ) 

+ -sin 2 0[x o 2 (— + —) + x, 2 (— + —)- ^ - ^] (30) 
2 R 0 R ] R l R 2 #j R 2 

+ -[y () 2 (- + -) + y, 2 (- + -)-^-^], 

2 R n R, R, R~ R, R~ 
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To find the equivalent C2 configuration for C3, the terms that contain the x 0 and x; 
coordinates in Eq. (30) should equal the corresponding terms in Eq. (13). This is trivial by 
setting 

x'j=xjsm0 (31) 

while maintaining all the inter-plane distances. 

The grating transmission functions rescale according to the general rule expressed in Eq. 
(24). In this case, substituting Eq. (31) into Eq. (24) results in the rescaled grating 
transmission functions for the equivalent C2 configuration: 

G' m (x m ,y m ) = G m (x m /sm8,y m ). (32) 
The detected pattern rescales according to 

I 3 (,x 2 ,y 2 )~ I' 2 (x' 2 ,y' 2 ), (33) 

where I 3 is the detected intensity pattern of the C3 configuration, and I 2 ' the detected pattern 
of the equivalent CI configuration. Substituting Eq. (31) into Eq. (33) leads to 

I 3 (x 2 ,y 2 )~ I' 2 (x 2 sm9,y 2 ). (34) 

The relationships expressed in Eq. (32) and (34) together describes an equivalent cone- 
beam, normal incidence interferometer for a cone beam, oblique incidence interferometer. 
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Fig. 3. Transformation of an extended source in the XZ plane as part of the process to map the 
oblique incidence configuration (C3) to an equivalent normal incidence configuration. With an 
extended source, the para-axial distance between any point on the source and the GO grating 
must be maintained between the two configurations. This condition is met by the shear 
transformation of the source in the pitch (XZ) plane. 

We now consider an extended area source such as an x-ray tube. Such a source can be 
viewed as an integral of mutually incoherent point sources. Each point on the extended source 
emits a cone beam, and the image intensity on the screen is the sum of the image intensities 
from all points. For simplicity we set the beam axes of all points to be parallel to each other. 
We then need to take into account that the para-axial distance between points on the extended 
source and the GO grating plane, R 0 , will change when the oblique GO grating is replaced by 
an equivalent vertical GO grating, except for points on the center horizontal line of the source. 
To maintain equality of the phase factors between the two configurations for all points on the 
source, it is necessary to eliminate the RO changes. This can be accomplished by moving each 
point on the source along its beam axis to offset the change. Therefore the area source 
undergoes a shear transformation in the XZ plane between the two equivalent configurations. 
This is illustrated in Fig. 3. 

5. From cone beam oblique incidence parallel planes to non-parallel planes 

The non-parallel plane C4 configuration is defined by 

C4 : (finite /? 0 ,(9 0 =9,9, =9 + S i ,9 2 =9 + S 2 ,G m ). (35) 
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The phase factors of C4 and C3 are given by Eq. (8) and Eq. (30), respectively. To find 
the equivalent C3 to C4, the phase factors must differ by no more than a constant, meaning 
that the corresponding terms containing the coordinates on the GO and Gl planes should be 
equal. Again the X coordinate terms will be considered first since they involve the incident 
angles. By choice we will keep the GO and Gl transmission functions unchanged, so that the 
equivalence conditions are 

x' 0 =x 0 , (36) 
x\ = x p (37) 

I » ( 7V ) = X ° 2( FV ) ' (38) 

*'o*'i -^(i + ^cotfl), (39) 



x 'i 2 ( i + IF* = x ' 2 + + 2<5 > cot ^ (40) 

"l K 2 K \ K 2 

^ = ^[l + (S l+ S 2 )cot0l (41) 
K 2 R 2 

Substituting Eq. (36) and (37) into Eq. (39) yields to the first order of S l 

R\ aR\(l-6 1 cot0). (42) 
Substituting Eq. (36) and (42) into Eq. (38) yields 

R' o *R o (l + ^-S l cot0). (43) 

Substituting Eq. (37) and (42) into Eq. (40) yields 

R\«R 2 [\-{2 + -^)5, cot 0] . (44) 
Finally, substituting Eq. (37) and (44) into Eq. (41) yields 

x\=x 2 [\-{\+ -^-)5 i cot 0 + 8 2 cot 0] . (45) 
Next let us consider the scaling of the Y coordinates. The requirements are 

y '' ( ^ + ^ ) = y ° 2( T + ^ (46) 

K 0 K 1 "() "l 



y o y i _ y 0 y t 



R\ R, 



y t y 2 _ yji 



R\ R 2 



(47) 
(48) 



y \ 2 (—+—) = y, 2 (— + — ), (49) 

K I K 2 K l K 2 
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y i y 2 _ ^ 



(50) 



Substituting Eq. (42) and (43) into Eq. (46) yields 

y'o = Jo- (51) 
Substituting Eq. (42) and (51) into Eq. (47) yields 

y\ = yi (l-S lC ot0). (52) 

Equations (42), (44) and (52) together also guarantee the condition of Eq. (49). Lastly, 
substituting Eq. (44), (52) into Eq. (50) yields the scaling of the y 2 coordinate: 

y' 2 = y 2 [l-(l + ^)S lC ot0l (53) 
K i 

We see that slightly tilting the Gl and detector planes relative to the GO plane is 
equivalent to changing the inter-plane distances according to Eq. (42) to (44). There is also an 
accompanying rescaling of the grating transmission functions according to the transformation 
of the (xj, yj) coordinates. Specifically, the grating transmission functions of the equivalent C3 
configuration should satisfy the general rule of Eq. (24), leading to 

G 'o( x (p3' () ) = G () (x 0 ,y () ), (54) 

and 

G \ (jq , y, ) = G, (jq , (1 + 8, cot 0) y, ). (55) 

The scaling of the coordinates on the detector plane given by Eq. (45) and (53) means that 
the detected intensity image of the non-parallel plane (C4) configuration is related to that of 
the parallel plane (C3) configuration by the relationship 

7 4 (x, , y 2 ) ~ / ' 3 ([1 - (1 + —)S l cot 0 + S 2 cot 0]x 2 , [1 - (1 + — )<?, cot 0] y 2 ), (56) 
R i R i 

where I 4 is the detected intensity pattern of the C4 configuration, and 7 3 ' that of the equivalent 
C3 configuration. Equations (42)-(44) and (54-56) together relate a cone beam, slightly non- 
parallel plane (C4) configuration to its equivalent cone beam, parallel plane (C3) 
configuration. 

Similar to section 4, we can extend the above concept to an area source, which is an 
assembly of mutually incoherent point sources. In the scaling relationships between the non- 
parallel plane (C4) and the equivalent parallel plane (C3) configurations, we note that only the 
rescaling of the source-to-GO distance R 0 is dependent on R 0 itself. Since different points on 
the extended source have different R 0 's, each source point needs to be moved para-axially by 
the appropriate amount to satisfy the rescaling relationship, resulting in a shear transformation 
of the source in the pitch plane, similar to the one illustrated in Fig. 3. The rescaling of R u R 2 
and the image on the screen are independent of R 0 and therefore remain the same as those of a 
point source. 

6. Summary of theoretical results 

All of the above results are summarized in Table 1 for convenience of reference. The 
following discussion on the effects of non-parallel grating planes will be based on the 
relationship between the C4 and the equivalent C3 configuration. 
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Table 1. Conversion Between Equivalent Interferometer Configurations 



Actual configuration 


Equivalent, simpler 
configuration 


Image on detector plane 


C2: cone beam normal incidence 


CI: Parallel beam normal 
incidence 




Ro 


CO 




R, 


R,(l +R,/Ro) 








ir(x 2 (Ro + Ri)/(Ro + Ri + 
R 2 ), y 2 (Ro + Ri)/(Ro + 
Ki -t- k 2 ;; 


R 2 


R 2 (Ro + R,)/(R 0 + Ri + R2) 


9„ = (90°, 90°, 90°) 


9„ = (90°, 90°, 90°) 




Go(x 0 , y 0 ) 


G 0 (xo/(l + R,/R„), y 0 /(l + 

t\.\/tS.())) 




Gi(x,,y,) 


Gi(x,,y,) 




C3: Cone beam oblique 
incidence 


C2: Cone beam normal 
incidence 




Rj 


R.i 


I 2 '(x 2 sin9, y 2 ) 


e„ = (e, e, 6) 


9„ = (90°, 90°, 90°) 




G m (x m , y m ) 


G m (x m /sm0, y 0 ) 




C4: Cone beam oblique 

!n/^!/~I/?vi/^/? M/iH tin iv i ]] pi nlnviPv 

If^ltiCfiLC, fllJlt fSLifLlllCl j/lLltlLO 


C3: Cone beam oblique 

inf^i/'lPn/'p n/i t'/ills> 1 til/i n o c 

It IL, lU,CflUC, fSLllLllLCL [SlUtLCd 




Ro 


R 0 [l + (Ro/Ri)5,COt9] 




R, 


R,(l-5iC0t9) 








I 3 '(x 2 {l-[(1+R 2 /Ri)5i- 
5 2 ]cot9},y 2 {l-(l + 
R 2 /Ri)5,cot9}) 


R 2 


R 2 [l-(2 + R 2 /Ri)6iC0t9] 


e n = (e, e + 5,, e + s 2 ) 


9„ = (9, 9, 9) 




Go(x 0 , y 0 ) 


Go(x 0 , y 0 ) 




Gi(x,,y,) 


Gi(xi, y,(l + 8iC0t9)) 





7. The effect of and correction for slightly non-parallel grating planes 

A useful insight from the theoretical modeling above is the influence of small pitch angle 
differences between the three principal planes. Based on the equivalence between the C3 and 
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C4 configurations (Table 1), the effect of tilting the detector plane is simply a rescaling of the 
image, and will therefore not be considered further. The effect of non-parallel grating planes 
is proportional to the cotangent of the incident angle. Thus it diminishes as the incident angle 
approaches 90°, and increases with decreasing incident angles. Therefore, non-parallelism 
deserves attention with oblique incidence gratings [11], and particularly in grazing angle 
reflective interferometers where the incident angle is necessarily smaller than the critical 
angle of total reflection on the grating surfaces [12-17]. In the following we will look at both 
the transmission and reflection types of interferometers. 

In the cone beam oblique incidence configuration C3, the leading-order effect is 
equivalent to changing of the various distances while maintaining parallel planes, as 
expressed in Eqs. (42)-(44). Those results imply that slightly non-parallel grating planes can 
be compensated for by adjusting the distances between the principal planes. Although there is 
also an accompanying rescaling of the Gl grating transmission function in the Y direction as 
expressed in Eq. (55), we will confine the discussion to linear gratings in the X direction (see 
Fig. 2), where rescaling in the direction out of the pitch plane does not change the detected 
image. An apparent question is why we prefer distance adjustments over angle adjustments? 
The reason, as detailed below, is that at small incident angles, a sub-milliradian or arc minute 
pitch angle difference translates to millimeter or centimeter positional adjustments. Looking 
at the required precisions of the two types of adjustments, it is apparent that the positional 
adjustments are easier to make than the corresponding angular corrections, particularly with 
compact stages. Thus, this type of compensation has practical advantages. 



6i Gl Detector 




Source 



Fig. 4. Compensating for a small difference eh between the pitch angles of GO and Gl gratings 
by moving the detector plane along the beam axis over a distance of AR 2 . The quantitative 
relationship between the pitch angle difference and the compensating movement of the detector 
plane is given in Eq. (66). 

In order to obtain the positional adjustments that compensate for a given pitch angle 
difference between the grating planes, we focus on the relationship between the C4 and 
equivalent C3 configurations, which is summarized in Table 1. The first step is to find the 
Talbot self-imaging conditions for the C3 configuration. The next step is to see how tilting the 
Gl plane disrupts the self-imaging conditions, and lastly the appropriate positional 
adjustments to restore the condition. 

The Talbot self-imaging conditions for the C3 configuration can be derived from that of 
the C2 normal incidence configuration, given by Donath and co-authors [19] 



R,R 2 1 (R, + R 2 ) = aP, 2 IX, 



(57) 



(R+R 2 )/R 2 =bP 0 /P, 



(58) 



Where a and b are the appropriate constants for the type of gratings used in the interferometer 
and a chosen Talbot order, P 0 and P 1 are the periods of the GO and Gl gratings, and "k is the 
wave length. Based on the results summarized in Table 1, the equivalent C2 configuration of 
the C3 involves only a rescaling of the grating transmission functions, which means that the 
grating periods in the equivalent C2 configuration are 



P\ =P sin0. 



(59) 
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Equation (59) reflects the intuitive idea that the effective grating period under oblique 
incidence is the projection of the actual grating period along the beam, and can thus be greatly 
reduced when the incident angle is small [12-17,20]. Substituting Eq. (59) into the self- 
imaging conditions of Eqs .(57) and (58) then yields the self-imaging conditions for the 
oblique incidence, parallel plane geometry: 

R l R 2 l(R x +R 2 )lsm 2 9 = aP 2 I X, (60) 

(R l + R 2 ) I R 2 = bP 0 I P v (61) 
Next we note that according to the relationship between the C3 and C4 configurations 
expressed in Eq. (42) to (44), only the tilt of the Gl plane and not the tilt of the detector 
screen changes the Rfs, of the equivalent C3 configuration, and thereby influence the fringe 
visibility. To quantify deviations of the self-imaging conditions caused by the Gl tilt, the 
equivalent R ) in Eq. (42) to (44) are substituted into the left-hand side of Eq. (60) and (61) to 
yield 

R[R' 2 l{R[ + R' 2 )l sin 2 9 = (1 - 2S 1 cot 9)aP 2 1 A, (62) 

(Rl + R' 2 )I R 2 =(l + S 1 cot 9)bP 0 1 P 1 . (63) 
Of the above two deviations, the first has less impact on the fringe visibility for two 
reasons: firstly the dependence of the fringe visibility on the parameter a in Eq. (60) is cos(Tia) 
[3], which is weak around the a value of maximal fringe visibility; more importantly, for a 
polychromatic source there is a continuous spectrum of Talbot distances, which substantially 
relaxes the condition of Eq. (60). In fact, with the push towards smaller grating periods and 
higher sensitivities, the Talbot-Lau type of interferometers evolve into Mach-Zehnder type of 
interferometers [10,17,21] that work outside the Talbot self-imaging regime. Therefore, we 
will focus on meeting the second condition expressed in Eq. (61), which is independent of 
either wavelength or incident angle. 

Having only one condition to satisfy also means that only one of the distances needs to be 
adjusted. The preferred choice is to move the detector and leave the gratings unchanged. 
Doing so avoids disturbing the grating pitch angles. This method is illustrated in Fig. 4. The 
pitch angle of the detector plane may be altered during its position adjustments, but this does 
not influence the self-imaging conditions expressed in Eqs. (62) and (63), and therefore does 
not affect the performance of the interferometer. 

To determine the correct amount of R 2 adjustment in Fig. 4, denote the adjusted distance 

as 

R 2A =(l + A)R 2 . (64) 
By substituting Eq. (64) into Eq. (44) to obtain the adjusted R 2 of the equivalent C3 
configuration, and then into the left hand side of the self-imaging condition in Eq. (63), 
retaining the leading order terms in the tilt angle d\, we have the adjusted form of the second 
condition: 

{R[ + R' 2A )I R' 2A =[1-ARJ (R, + R 2 W + 8, cot 9)bP a I P v (65) 

By setting the right hand side of Eq. (65) to the same as Eq. (61), we obtain the 
appropriate adjustment AR 2 that restores the self-imaging condition: 

AR 2 = R 2 8 x cot 9(R l +R 2 )/R 1 . (66) 
As an example of R 2 distance compensation for a pitch angle difference between the Gl 
and GO gratings, we consider a symmetric interferometer of Ri = R 2 = 22.5 cm. The amount of 
R 2 adjustment as a function of the tilt angle is given by Eq. (66) and graphed for several 
incident angles 6 in Fig. 5. In this example, a 5 arc minute or 2.9 milliradian Gl pitch error is 
compensated for by R 2 adjustments of 0.38 to 1.1 millimeters depending on incident angle. 
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Fig. 5. The relationship between the pitch angle error of the Gl grating (relative to the GO 
grating) and the translational movement of the detector plane that compensates for that error. 
The three curves are for three different x-ray beam incident angles of 30°, 45° and 60°. Other 
parameters of the interferometer are given in the text. 



Generally compact goniometers are not graduated below 5 arc minutes, while micrometer 
driven linear stages have graduations at 0.01 millimeters. So it is practically easier to adjust 
the detector position than the grating pitch angle. The example also shows the trend towards 
larger corrections as the incident angle becomes smaller, or more oblique. 

Lastly we consider applying the above compensation method to the case of an extended 
source, which is viewed an assembly of mutually incoherent point sources. Instead of setting 
the beam axes for all point sources to be parallel as in sections 4 and 5, here we define the 
beam axis from a source point as the straight line connecting the point to the pivot axis of the 
Gl grating and parallel to the pitch plane (Fig. 4). Now each point on the source has its own 
incident angle. Because the lateral adjustment of the screen (Eq. (66)) is dependent on the 
cotangent of the incident angle, the compensation method will only work properly if the range 
of incident angles subtended by the source is small relatively to the mean incident angle. In 
geometric terms the condition is that S/(R 0 + Ri) « 8, where S is the vertical size of the 
extended source. 

8. Application to x-ray reflective interferometers 

X-ray reflective interferometers (Fig. 6) have small incident angles by nature of the total 
reflection requirement. A reflective intensity grating is a surface pattern of periodic, 
alternating stripes of high and low reflectivity coatings [12-16] such as gold stripes on glass. 
A reflective phase grating has periodic elevated stripes to cause oscillating phase delays [17] 
on the reflected wavefront. X-ray interferometers of this type operate on the premise that the 
surface materials have refractive indices that are less than 1.0, such that total internal 
reflection occurs in air which has a refractive index of 1 .0. 

Reflective interferometers can be modeled by their equivalent transmission configurations 
[17] as illustrated in Fig. 6. Here the waves and optical components down-beam from the GO 
grating are represented by their mirror images with respect to the GO surface, and the waves 
and components down-beam from the Gl grating are further reflected with respect to the Gl 
surface. The appropriate compensation for a slight tilt of the Gl grating has been calculated 
for the transmission equivalent in the previous section, and so, by the mirror reflection 
operations, the same information for the reflective configuration can be obtained. 
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><T2, ZT2I 



(x't2,z't 2 ) 



Fig. 6. Illustration of a grazing angle reflective grating interferometer in black line drawings, 
and its equivalent transmission configuration (ETC) in red line drawings. In the ETC the 
detector is represented by its reflection with respect to the Gl grating surface, and both the 
detector and the Gl grating further reflected with respect to the GO surface. In the real 
reflection configuration, the intersections of the beam axis with the Gl grating and the detector 
plane are (xi, z\) and (x 2 , Zi), while the equivalent points in the ETC are (x T i, Zti) and (jt T 2, Zti)- 
To find the new detector position (ac \ z'%) that compensates for a small pitch angle difference 
5\ between the Gl and GO gratings, we first find the new detector position (x't2, z'n) in the 
ETC according to Eq. (66), then reflect it with respect to the GO surface and the tilted Gl 
surface. 

To carry out this procedure, the pitch plane of the interferometer is set as the XZ plane, 
with the X axis parallel to the GO grating plane (Fig. 6). The center of the x-ray beam 
footprint on the GO surface is set as the origin of the XZ coordinates, and the beam centers on 
the Gl and detector planes have the coordinates (x u z\) and (x 2 , z 2 ). Their corresponding 
points in the equivalent transmission geometry are labeled by (x T1 , z T1 ) and (x T2 , z-n)- When 
the Gl plane is tilted by a small angle S t , the compensating position of the detector in the 
transmission geometry is labeled by (x' T2 , z't2), and the corresponding detector position in the 
reflective configuration is (x \, z ' 2 ). Other parameters of the interferometer are the distances 
between GO and Gl of R u and between Gl and the detector of R 2 , and the grazing incident 
angle 9. Referring to Fig. 6, in the absence of the Gl tilt the various coordinates are 



When the Gl plane is tilted, the appropriate adjustment of the R 2 distance to compensate 
for the tilt in the transmission geometry is given by Eq. (66). Substituting into Eq. (67) yields 
the adjusted coordinates (x' T2 , z'ri)'- 



where AR 2 is given by Eq. (66). By reflecting with respect to the tilted Gl surface and then 
the GO surface, the adjusted position of the detector in the actual reflection geometry is given 



(x 1 ,z 1 ) = (R 1 cos 0, R t sin 0), 

(x 2 ,z 2 ) = (R 1 cos 0 + R 2 cos 0, R l sin 0 — R 2 sin 0), 

(x Tl ,z n ) = (Ri cos 0, sin 0), 

(x T2 ,z T2 ) = (Ri cos0 + R n cos 0, sin0-R 1 sin 0). 



(67) 



(x' T2 ,z' T2 ) = (Ri cos 0 + R 2 cos 0 + AR 2 cos 0, 
- R t sin 0 - R z sin 0 - AR 2 sin 0), 



(68) 



by 
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(x' 2 ,z' 2 ) = (x 2 +Ax 2 ,z 2 +Az 2 ), 
Ax 2 = S l R 2 (l + R 2 / RJcotO, 
Az 2 = -3^(3 + R 2 /R 1 ) cos 0. 



(69) 



To illustrate the detector adjustment in practical settings, we consider x-ray reflective 
gratings that use a high density metal coating such as gold or tantalum to form the grating 
pattern [13,16]. The grazing incidence angles are in the range of milliradians. In an example 
of Ri = R 2 = 22.5 cm, Fig. 7 plots the lateral and vertical compensations of the detector 
position for a range of Gl pitch angle errors at several incident angles. The graph shows that a 
100 micro radian (21 arc seconds) pitch angle error translates to a lateral adjustment of the 
detector by 15 to 45 mm, depending on the incidence angle. In contrast, the vertical 
adjustment is only 0.09 mm and it is independent of the incident angle. In practical terms, it is 
far less demanding to make millimeter adjustments of the detector position than arc second 
adjustments of the grating planes if we use off-the-shelf commercial compact stages. 
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Fig. 7. The relationship between the pitch angle error of the Gl grating (relative to the GO 
grating) and the movement of the detector plane that compensates for that error in a grazing 
incidence reflective interferometer, (a) The required axial movement of the detector as a 
function of the pitch angle error for three different grazing incident angles of the x-ray beam, 
(b) The required vertical movement of the detector as a function of the pitch angle error does 
not change with the incident angle. Other parameters of the reflective interferometer are given 
in the text. 



9. Experimental demonstration in a hard x-ray reflective grating interferometer 

We built a hard x-ray reflective grating interferometer of the type illustrated in Fig. 6. The 
interferometer consisted of a compact tungsten anode x-ray tube operating at 80 kVp/lmA, a 
GO reflective intensity grating which was a 5.0 cm diameter flat glass surface coated with 30 
nm thick gold stripes at a period of 2.1 mm, a Gl reflective phase grating of gold coated flat 
glass surface with elevated stripes (5 nm elevation) also at 2.1 mm period [17], and a GdOS 
phosphor screen (Beam Imaging Solutions Inc) serving as the detector plane. The grazing 
angle was set to 2.0 milliradians. The grating stripes were all in the Y direction perpendicular 
to the pitch plane (Fig. 6). When observed along the beam at the grazing angle, the effective 
grating period was 4.2 |im. In terms of the contrast of x-ray reflection between gold-coated 
and bare glass surfaces, the former was approximately 8 times the latter when measured by 
the light intensity from a phosphor screen. The distance parameters were set at R 0 = 25cm, 
= R 2 = 22.5 cm. The gratings and the phosphor screen were all mounted on compact 
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goniometers with manual micrometer drives. The smallest graduation on the micrometers 
corresponded to an angular increment of 0.4 milliradians. 

No filter was used for the x-ray tube, and thus the beam had a broad spectrum. The Talbot 
self-imaging condition expressed in Eq. (60) was satisfied at two discrete energies in the 
spectrum, namely at 61 keV when the effective propagation distance was l/16th the Talbot 
distance, and at 21 keV at 3/16th Talbot distance. Because the phosphor screen was 
illuminated at grazing incidence, the fringes were expanded and could be resolved with a 
standard Nikon D3X camera. Top-down images of the phosphor screen were acquired at 60 
sec exposure and matrix size of 4256 x 2832. The pixel size on the plane of the screen was 
determined by two markers on the screen of known spacing to be 15.7 pm. Ideally the 
phosphor screen should be parallel to the gratings, however we found that the x-ray flux was 
too dispersed on the screen to give sufficient phosphorescent light in 60 sec exposure. Instead 
we increased the pitch angle of the screen to 5.6 milliradians. 
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(c) 

Fig. 8. X-ray reflective interferometer demonstration of compensating for a small pitch angle 
difference between the GO and Gl gratings by moving the detector phosphor screen, (a) In the 
initial setup, the image on the detector screen consists of the primary interferences fringes on 
the left, and large-period secondary fringes (alignment aids) on the right margin. The intensity 
profile of a center cross-section of the primary fringes is shown. The fringe visibility as defined 
in the text is 12.3%. The scale bar indicates 20 um in the vertical direction perpendicular to the 
x-ray beam, (b) After making a small change of the Gl pitch angle, the primary fringes are 
obliterated, (c) Moving the detector screen by 4.2 cm along the x-ray beam axis recovers the 
primary fringes, and the fringe visibility returns to 12.9%. The movement distance indicates 
that the the pitch angle change we introduced was 0.19 milliradians (39 arc seconds). 
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For initial alignment of the gratings, we used a simplified version of a method described 
by Liu and co-authors [15]. We coated on both reflective gratings a secondary grating pattern 
of a very large period alongside the primary grating pattern. On the GO grating this was a 20 
mm period stripe pattern, and on the Gl grating it was a 10 mm period pattern. When setting 
up the interferometer, the secondary gratings provided large and easily visible fringes on the 
phosphor screen. The target of the setup procedure was to bring the period of the large fringes 
to 40 pm in the vertical (Z) direction, which was an indication that the grazing angles were 
approaching 2.0 milliradians. After this step was accomplished, we then optimized the fringe 
visibility of the primary gratings by moving the phosphor screen in the longitudinal (X) 
direction. Figure 8a shows the image on the phosphor screen after the initial setup. The 
visibility of the primary interference fringes, defined as the ratio of (fringe peak - fringe 
trough)/(2*mean intensity), was 12.3%. The vertical width of the reflected beam was 
determined by the width of the Gl grating, which was 50 mm in diameter, or 100 |im in the 
vertical direction under the 2.0 milliradian incidence angle. This is in agreement with the 
width of the intensity profiles in Fig. 8. 

Next we made a small change of the pitch of the Gl grating to observe the effect on the 
interference fringes. We adjusted the goniometer of Gl by a fraction of the smallest 
graduation of 0.4 milliradians, but not known precisely. The fringe visibility dropped from the 
initial 12.3% to 1.3%, which is illustrated in Fig. 8b. We then attempted the distance 
compensation method by moving the phosphor screen away from the Gl grating, until the 
fringe visibility was recovered to 12.9% (Fig. 8c). The optimal distance to move the screen by 
was 4.2 cm. Based on the calculation in Eq. (69), the Gl pitch change which we introduced 
was 0.19 milliradians. 

10. Discussion 

We showed from first-principle wave propagation modeling that oblique incidence, cone 
beam illuminated interferometers with either perfectly parallel planes or slightly non-parallel 
planes have their equivalent parallel beam, normal incidence configurations if the Fresnel 
approximation is valid. The theory also shows that oblique incidence interferometers become 
increasingly sensitive to non-parallelism between the grating planes as the incident angle 
becomes smaller. The effect is particularly noticeable in x-ray reflective grating 
interferometers, where the incident angles are on the order of milliradians. Aligning the 
grating planes by directly adjusting their pitch angles can be a demanding task with compact 
stages. Fortunately, the equivalence relationships imply that the pitch angle errors can be 
compensated for by changing the detector position. In reflective interferometers, arc second 
pitch errors translate to millimeter positional adjustments, which can be realized easily. 

Although an ideal monochromatic Talbot Lau interferometer provides 100% fringe 
visibility, our observed values were 12 - 13%. The primary factor that lowered the fringe 
visibility was the broad spectrum of the source, which had two consequences. The first is that 
the Talbot self-imaging condition was only met at two discreet energies in the spectrum while 
the rest of the spectrum did not satisfy the Talbot condition. The second is that the phase 
modulation from the Gl phase grating was wavelength dependent [17], and therefore was not 
the targeted value of 180° over the whole spectrum. In contrast, a grazing angle far- field 
Mach Zehnder interferometer can be truly achromatic as shown in ref [17], and achieve 
substantially higher fringe visibility. Another factor that likely affected fringe visibility was 
the point response function of the phosphor screen. It was a 15 pm layer of GdOS powder of 
5 pm nominal grain size on a glass surface. The associated surface roughness was not 
negligible relative to the 4.2 |im fringe period, and may likely have blurred the fringes when 
used at grazing incidence. 

An additional observation of the results in Fig. 8 is that the fringe visibility tends to be 
decrease toward the top and bottom edges of the reflection band on the screen. This can be 
explained by the relatively large deviation of the phosphor screen from being parallel to the 
grating planes. In order to have sufficient x-ray flux density on the screen, the pitch angle of 
the screen was set at 5.6 milliradians. The corresponding pitch angle error S 2 (the deviation 
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from being parallel to the grating planes) was 3.6 milliradians. Although in the leading order 
approximation S 2 does not affect fringe visibility, this error was nearly twice the incident 
angle. Consequently, when the condition of Eq. (61) was met at the center of the screen, it 
was no longer met towards top and bottom edges of the reflection band. This notion was 
supported by our observation that when moving the screen towards the optimal position, the 
fringes first appeared at one edge of the reflection band and then moved to the center. Overall, 
the use of a phosphor screen at grazing incidence has the advantage of direct visualization of 
the interference fringes, but with inherent limitations associated with dispersed x-ray flux and 
surface roughness. An alternative, indirect imaging of the fringes is to replace the screen with 
another reflective intensity grating as the "analyzer" [4-6,16], followed by a vertical phosphor 
screen to detect the x-rays reflected from the analyzer grating. 

Because the Fresnel approximation is the starting point of this work, it is worth noting that 
for oblique incidence geometry, the Fresnel approximation is valid when the third order term 
in the Taylor expansion of the phase factor satisfies 

^cot<9«l, (71) 

where X is the x-ray wavelength, P is the typical grating period, R is the typical distance 
between the principal planes of the interferometer, and 0 is the incident angle. The third order 
term vanishes for normal incidence, but can become significant for small incident angles. The 
associated aberration effect for Talbot fringes has been analyzed by Testorf and co-authors 
[22]. For our reflective x-ray interferometer it was approximately 0.004 and negligible. 
Ultimately, when the grating period becomes sufficiently small, the higher order terms will be 
significant and we transition from a near-field Talbot Lau interferometer to a far field 
situation. Nevertheless, in such cases coherent interference effects can still be obtained with 
either normal or oblique incidence using polychromatic incoherent sources, as long as certain 
symmetry conditions are met. This idea has been theoretically analyzed and experimentally 
demonstrated in matter wave and light wave interferometers [1,2,10,17]. 
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